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Abstract 

Kusuoka and Zhou have defined the Laplacian on the Sierpinski carpet 
using average values of functions on small cells and the graph structure of 
cell adjacency. We have implemented an algorithm that uses their method 
to approximate solutions to boundary value problems. As a result we have 
a wealth of data concerning harmonic functions with prescribed bound- 
ary values, and eigenfunctions of the Laplacian with either Neumann or 
Dirichlet boundary conditions. We will present some of this data and 
discuss some ideas for defining normal derivatives on the boundary of the 
carpet. 
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1 Introduction 

The Sierpinski carpet, SC, is the self similar fractal in the 
the identity 
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SC=\jF i (SC) 

where {Fi} are contraction mappings defined by 

Fi(x) = l/3(> - Qi) 
where {qi} are the eight points shown in Figure [Ll] 
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Figure 1.1: The level 1 graph approximation of SC with the eight cells labeled. 

By convention, the boundary of SC, 9SC, is defined to be the boundary of 
the unit square containing all of SC. 

We follow the usual convention for iterated function systems like {Fi}, let 
uj = (cji, ...,u; m ) denote a word of length m whose each ujj takes on values in 
{0,1,2,3,4,5,6,7} and F u = F Ul o F^ o • • • o F Um . We call F W (SC) an m-cell, 
and let SC m denote the graph whose vertices are the m-cells and with edges 
connecting two m-cells that intersect along an edge. 

Originally there were two different constructions of a self-similar energy and 
Laplacian on SC with the full D4 group of symmetries, one due to Barlow and 
Bass [BB1993 and one due to Kusuoka and Zhou |KZj . It was recently shown 
that the two constructions are equivalent [BBKT] , and in fact the energy and 
Laplacians are unique up to a constant multiple. In particular this implies that 
the graph Laplacians on SC m , when appropriately normalized, converge to the 
Laplacian on SC m (in [KZ this was only proved for a subsequence of ra's). In 
this paper we systematically use this observation to approximate many of the 
basic functions on SC: eigenfunctions of the Laplacian (and their eigenvalues) 



plane determined by 

(1.1) 

(1.2) 
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under various boundary conditions, harmonic functions with prescribed bound- 
ary conditions, Poisson kernels, heat kernels, and Dirichlet kernels. (A different 
approximation method was used in [BHS , but it does not appear to be as ver- 
satile.) In addition we study analogous problems on fractafolds and covering 
spaces of SC. For example, one of our most striking experimental observations 
is that the fractafolds obtained from SC by gluing the boundaries according 
to the torus, Klein bottle, or projective plane model, have eigenvalue counting 
functions that differ from each other by uniformly bounded amounts. This is a 
stronger statement than is true for the standard Laplacian on the torus, Klein 
bottle, and projective plane. 

We note that SC is quite different from the PCF fractals ( |Ki[ S3 ) such as 
the Sierpinski gasket SG, where there is a simpler construction of a Laplacian. 
In particular, for SG there is a theory of spectral decimation ( |FSj ) that allows 
the exact computation of eigenfunctions of the Laplacian. In absence of precise 
theoretical tools on SC, the experimental approach becomes even more impor- 
tant. We will discuss the implications of our numerical results in the concluding 
section of this paper. 

We present a brief discussion of the method of averages in section 2. In 
section 3 we present our results on harmonic functions, including a discussion of 
the Poisson kernel and the effective resistance metric. In section 4 we discuss the 
spectrum of the Laplacian under Dirichlet and Neumann boundary conditions, 
and also for fractafolds [Si] obtained by identifying the boundary as in the 
construction of the torus, Klein bottle, and projective plane from the square. 
We discuss many properties of the eigenvalues, and also present some graphs 
of the eigenfunctions. The actual eigenfunctions are crucial to the subsequent 
computations in this paper. In section 5 we discuss the heat kernel, both the 
on-diagonal and off-diagonal behaviors. In section 6 we examine the decay rates 
at the boundary of some of the functions we have constructed. Section 7 is 
devoted to the spectrum of two examples of covering spaces with a Z-action, in 
the spirit of results for covering spaces of Julia sets in |ADSj . Section 8 contains 
a concluding discussion. 

One of the motivating problems for our work was the challenge of trying to 
define a normal derivative on the boundary of SC for functions in the domain of 
the Laplacian in such a way that a Gauss-Green formula is valid. (Indeed this is 
one of the basic tools for studying Laplacians on PCF fractals.) In section 6 we 
discuss our experimental observations that indicate why this is a very challenging 
problem. So, although our work does not resolve this problem, it may provide 
guidance for future work. We present a small sampling of our numerical results 
here. More data, and the programs used to generate them, may be found on 
the website |http: //www. math. Cornell . edu/~reu/ sierpin ski- car pet/. 

2 Method of Averages 

The method of averages was introduced by Kusuoka and Zhou in |KZj for func- 
tions on SC. Instead of function values defined at vertices in fractal approxima- 
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tion graphs, we consider average values on entire m-cells. Let {A^} be the 
collection of all m-cells. Then we define f(A^) = 8 m f Am) f dfi. 

J 3 

We reference a specific m-cell by the cell address which is the sequence of 
integers corresponding to the contraction maps that map the unit square to the 
particular cell. For notation, we list the cell address in single quotes and in the 
opposite order that it would be read as function composition. For example, if / 
is the unit square, then the 3-cell given by F o F 3 o F 2 {I) is '230'. 

We write ~ if the two cells have a boundary edge in common. 

We say that A^ and A^ are neighbors. It can be verified that any cell in 
SC m will have either 2, 3, or 4 neighbors. There exist algorithms that give 
the number and specific cell addresses of the neighbors of any m-cell and are 
available on the website. 

Using this method of averages we can now define discrete approximations 
using graph energy and graph Laplacian on level m: 

E m u)= E (/(^ m) )-/(4 m) )) 2 . (2-i) 

A m f(A™) = E (/(4 m) )-/(4 m) ))' ( 2 -2) 

X (»W(») 

To compare these approximations on different levels we need to introduce 
renormalization factors. For the energy there is a factor r such that 

£ m (f) = r- m E m (f) (2.3) 

and 

£(f) = lim E m (f) (2.4) 

constructs a local, regular Dirichlet form (energy) on SC. Experimentally, r _1 ~ 
1.25. Then 

Af = lim (Sr-^A^f (2.5) 

m— ^oo 

is the Laplacian on SC. We denote by p the factor 8r _1 ~ 10.01. The energy 
and Laplacian are related by 

/ (Af)gdn = £(f,g) (2.6) 

JSC 

for g G dom£ with g vanishing on dSC (here £(f,g) is obtained from £ (/) by 
polarization), where \i is the symmetric self-similar measure on SC. The factor 
8 in ( |2.5| ) is of course the scaling factor for ja. 

We will be imposing various boundary conditions on SC. Since the boundary 
of SC is the boundary of the outer most unit square, we have to handle cells that 
intersect with the boundary of SC separately. To encode the behavior on the 
boundary to our method of averages on cells we introduce the notion of virtual 
cells. Every cell that borders the boundary has exactly one virtual cell on the 
exterior of SC except for the four corner cell which each have 2 virtual cells. 
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Figure 2.1: Every cell in SC will have exactly two to four neighbor cells. For 
example, on the left blowup of SC, the interior cell a has neighbors 6, c, and 
d. On the right hand picture, x is the top-left corner cell. It has four neighbor 
cells: two interior cells, y and z\ and two virtual cells x and x' . 



3 Harmonic Functions 

In this section we study harmonic functions on SC by approximating them 
by harmonic functions on SC m . At interior cells x, the harmonic condition is 
simply that h(x) is the average value of h(y) over the neighbors y (there may 
be two, three, or four neighbors). We would like to impose boundary conditions 
h\dsc = f where / is a function defined on <9SC. It would be clumsy to try 
to impose these conditions on cells in 9SC m , since these cells extend into the 
interior of SC. Instead we introduce "virtual cells" for each boundary cell x 
in SC m that reflect outside SC across the boundary line segment (for the four 
corner cells there are two such reflections, but we will not introduce any new 
notation to handle these cases). The function h will be extended to be defined 
on the reflected cells. We then require the harmonic condition 

(K x ) - K x )) = for each x e SC m (3.1) 



including virtual neighbors and the boundary conditions become 

— (h(x) + h(x*)) = Average value of / on x D x* (3.2) 



for every virtual cell x* . Note that together (3.1 ) and (3.2 ) give the same number 
of equations as unknowns (8 m +4-3 m ), and it is easy to see that there is a unique 
solution for any given /. If we wish to impose Neumann boundary conditions 



on all or portions of the boundary, we replace (3.2) by the even extension 



h(x*) = h(x) (3.3) 



for the appropriate cells. Note that this has the effect of imposing (3.1) without 
adding the virtual cell neighbor. 

For our first example, we consider the boundary data f(t) = s'mnkt along 
one edge of the boundary and zero on the other three edges. By taking infinite 
linear combinations of these functions and their three rotations we could obtain 
any harmonic function. In Figure [3TT] we display the graphs for k = 1,2, 3, 4, 5, 6 
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(computed at level m = 6). In Figure 3.2 we display the graphs of the restrictions 
of these functions to the line one third of the way to the opposite boundary edge. 



We do not know whether or not these harmonic functions have finite energy 



on SC. By examining the renormalization energies (Table 3.1) for levels m = 
1,2,3,4,5,6 we see that it is reasonable to conjecture that they do have finite 
energy. In Figure |3.3| we show the graph of £ q as a function of k together 



with a linear approximation. It is difficult to speculate on the dependence for 
large values of k. Note that under the assumption that these functions have 
finite energy, it would follow by standard Fourier series arguments that they are 
orthogonal in the energy inner product for a fixed boundary edge (but not so for 
rotations). This would allow us to compute energies for ^a^hk if hk denotes 
these harmonic functions for a fixed boundary edge. 



m 


k = l 


k = 2 


k = 3 


fc = 4 


k = 5 


k = 6 


1 


0.4012 


0.2487 


0.4688 


0.0622 


0.0160 





2 


1.2709 


2.2158 


2.5081 


2.4536 


2.1673 


1.7703 


3 


1.5559 


3.3898 


5.2038 


6.8267 


8.2068 


9.2851 


4 


1.6286 


3.6863 


5.9569 


8.2735 


10.6618 


13.0465 


5 


1.6482 


3.7626 


6.1398 


8.6161 


11.2343 


13.9264 


6 


1.6535 


3.7835 


6.1890 


8.7067 


11.3823 


14.1488 



Table 3.1: Renormalized energies £ m ( u k) of the harmonic functions given in 
Figure 13.1 



For our second example we compute approximations to the Poisson kernel 
P(x, t) for x GSC and t G <9SC. This is the harmonic function in x with boundary 
values equal to the delta measure at t. This has a singularity as x approaches 
t. The Poisson kernel allows us to construct any harmonic function from its 
boundary values via the Poisson integral 

h(x) = ( P(x,t)f(t)dt. (3.4) 

JdSC 

We refer to the point t as the stimulus point. We approximate the delta 
function in the obvious way, taking the right side of ( |3.2| ) to be 3 m on the cell 
for which t G x U x* and zero for all other boundary cells (in case t lies at the 
boundary between two cells x\ and x^ we use 1/2 • 3 m for each of them). Our 
intuition is that the SC Poisson integral should resemble the upper half-plane 
Poisson kernel 

1 V 

P((x,y),t) = -7 ^ k here x, t G R and y > (3.5) 

7T (x — ty + y z 

at all points except the corner points, where it should resemble the quadrant 
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Figure 3.1: Graphs of the harmonic functions on SC6 with f(t) — sin7rH defined 
along one edge for k = 1, 2, 3, 4, 5, 6 and zero on the other three edges. 
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Figure 3.2: Restrictions of the harmonic functions from Figure 3.1 along the 
line one third of the way from the nonzero boundary edge. 



(x > 0, y > 0) case where the analog of (3.4) is 

Kx,y)= l -Y^l ( x i- y ^Tr + (2 Xy y f±{s)ds (3 - 6) 

where the boundary values are 

h(x,0) = f + (x) and h(0,y) = f_(y). (3-7) 
The two basic questions we would like to discuss are the following: 

(i) What do the level sets of P look like (for fixed stimulus point)? Do they 
resemble circles tangent to the stimulus point? 

(ii) What is the rate of growth as x approaches the stimulus point? 

We display the graphs of approximations to P(-,t) for three choices of stimulus 
point in Figure [3~4| The choices are the midpoint of an edge, a point one third 
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Figure 3.3: Graph of £s(uk) with linear approximation. 



of the way along an edge (the intersection of two 1-cells), and a corner point. 



In Figure |cT5] we display some thickened approximations to level sets for these 
choices, giving a pictorial answer to (i). To study question (ii) we restrict the 
Poisson kernel to a straight line segment that approaches the stimulus point 
from inside SC. In the case of the junction point stimulus point there is a line 
segment perpendicular to the edge, while for the center and corner points, there 
is a "half-diagonal" segment joining a corner point to an opposite center point. 



In Figure 3.6 we display the log of the restriction of the Poisson kernel to the 
line segment together with a straight line fitted to the data in the relevant 
portion of the graph (not too close to the stimulus point where the accuracy of 
the approximation deteriorates, and not too far from the stimulus point). We 
measure the distance to the stimulus point in the Euclidean metric. We observe 
a power law blow-up rate, but the power is different for each of the points. 

For our third example we use harmonic functions to compute effective re- 
sistances. The effective resistance metric R(x,y) between points x, y G SC is 
defined by 

R(x J y)~ 1 = mm{£(u) : u(x) = and u(y) = 1}, (3.8) 

and indeed the minimum is attained by a function that is harmonic in the 
complement of the points x and y and satisfies Neumann boundary conditions. 
For our approximations we can use the same formula for x and y cells in SC m 
and in fact we can interpret the result as the effective resistance between these 
cells in SC. 



In Figure 3.7 we show the graph of R(-,y) where y is the center point of 



a boundary edge. In Figure 3.8 we show thickened level sets of this function. 



These are clearly different from the level sets of the Poisson kernel. In Section 
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Figure 3.5: Thickened level sets to P(-,t) for stimulus point on the center, 
junction point, and corner of an edge on the boundary. 
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Center Poisson Kerne! restricted to half-diagonal line containing stimulus. Decay rate=-1 .1 575 




Figure 3.6: Log graphs of the restriction of P(-,t) to line segments for stimulus 
point on the center, junction point, and corner of an edge on the boundary. The 
center and corner kernels are along the "half diagonal" while the junction kernel 
is along the full perpendicular line segment. 
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5 we will compare them to level sets of the heat kernel. In Figure [3^9] we graph 
the restriction to the half- diagonal line segment and estimate the decay rate as 
x->y. 



4 Spectrum 

In this section we discuss the eigenvalues and eigenfunctions of the Laplacian 
with either Dirichlet or Neumann boundary conditions, and later other types of 
periodic boundary conditions. We note that a different method to approximate 
the Neumann spectrum was used in |BHSj . The results obtained using the two 
methods are consistent. In our method we solve the eigenvalue problem on SC m 

- A m u(x) = Xu(x) (4.1) 

at all cells in SC m , computing the Laplacian taking into account virtual cells 
using even or odd extensions 

u(x*) = ±u(x) (4.2) 

on virtual cells for Neumann or Dirichlet boundary conditions. The D4 sym- 
metry group acts on SC and all approximations, so each eigenspace must split 
according to the irreducible representations of D4 (four 1-dimensional and one 
2-dimensional). Following the notation in [BHS we denote the 2-dimensional 
representation by 2, and the 1-dimensional representations by 1 ++, 1 + — ,1 — + 

and 1 . The first ± refers to the symmetry (+) or skew- symmetry (— ) with 

respect to the reflections about the diagonals of the square, and the second ± 
refers to the symmetry or skew-symmetry with respect to horizontal and ver- 
tical reflections of the square. This implies that there must be eigenspaces of 
multiplicity two. In principle there could be higher multiplicities due to coinci- 
dences, but we do not observe any in our data, and there is no reason to believe 
that they occur. 

As discussed in [BHS], there is a miniaturization principle: for each eigen- 
function u with eigenvalue A there exists an eigenfunction u' with eigenvalue pA, 
with u r composed of eight miniaturized versions of u in each 1-cell, depending 
on the associated representation. It is easy to recognize the miniaturized eigen- 



functions and to use the eigenvalues to estimate the value of p. In Figure 4.1 



we show graphs of typical eigenfunctions with Dirichlet boundary conditions, 



and in Figure 4.2 the same with Neumann boundary conditions. In Tables 4.1 
and 4.2 we show an initial portion of the Dirichlet and Neumann eigenvalues 
on levels m = 5 and 6, along with the representation types of the associated 
eigenfunctions. We observe the "quartet" structure of each spectrum: if we 
group together four consecutive eigenvalues then each quartet contains a multi- 
plicity two eigenspace corresponding to the 2-dimensional representation of D4, 
and two distinct 1-dimensional representation eigenspaces. (there appear to be 
some short term exceptions to this structure in the Neumann spectrum, but it 
is not clear if these represent a true deviation to the quartet structure, or are 
merely due to inaccuracy in the computation higher up in the spectrum.) 
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Figure 3.7: Graph of R(-,y) on SC4 where y is the center-most cell on the top 
boundary edge. 
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Table 4.1: The first 60 renormalized Dirichlet eigenvalues including the unrenor- 
malized approximations for levels m = 5 and 6 along with the D4 symmetry 
type of the associated eigenfunction 
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Table 4.2: The first 60 renormalized Neumann eigenvalues including the un- 
renormalized approximations for levels m = 5 and 6 along with the D4 symme- 
try type of the associated eigenfunction 
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R( .,y) restricted to ha If-d iago nal line containing stimulus. Decay rate= . 3 S 5 2 




-3.5 -3 -2.5 -2 -1,5 
log(distance from stimulus point) 



Figure 3.9: Restriction of R(-,y) along the half diagonal on SC4 where y is the 
center-most cell on the top boundary edge. 



We computed the L°° norm of the eigenfunctions tpj (normalized so \\(pj\\2 — 
1). Figure 4.3 shows the values plotted as a function of the eigenvalue. The 
data is consistent with the uniform boundedness of H^/Hoo as j ~^ 00 ? al- 
though we could not realistically expect to detect slow growth rates such as 
1 1 1 1 00 = 0(log Xj). Of course uniform boundedness hold for the ordinary Lapla- 
cian on the square, but is in sharp contrast with the power law growth on general 
Riemannian manifolds and on SC (localized eigenfunctions provide obvious ex- 
amples). 



In Figures 4.4 and |4.5|we graph the eigenvalue counting functions 



N(t) = : Xj < t} 

on a log- log scale, confirming the predicted slope of 

log£ 



logp 



0.9026. 



In Figures |4.6| and |4.7| we graph the Weyl ratios 

N(t) 



W(t) 



(4.3) 



(4.4) 



(4.5) 



It is believed that asymptotically W(t) should approach a multiplicatively 
periodic function (period p) bounded and bounded away from zero (the same 
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Figure 4.2: The first nine nonconstant Neumann eigenfunctions on SC5 
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Figure 4.3: The values of max \cpj \ plotted as a function of Xj for levels 3, 4, and 
5 respectively. The plot for level 5 is incomplete as we are only able to compute 
the first 1500 eigenvalues but it still appears to exhibit the same behavior. 



Dirichlet (level 5) eigenvalue counting function 




t 



Figure 4.4: Eigenvalue counting function for the level 5 Dirichlet eigenvalues on 
a log- log scale. 



Neumann (level 5) eigenvalue counting function 




t 



Figure 4.5: Eigenvalue counting function for the level 5 Neumann eigenvalues 
on a log- log scale. 
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Dirich let (level 5) Weyl Ratio 




10" 3 10" £ 10" 1 10° 10 1 



t 

Figure 4.6: Weyl ratio for the level 5 Dirichlet eigenvalues. 



Neumann (Ur-J<3l 5) Weyl Ratio 



10 " - 




10" 3 10" 2 10" 1 10 a 



Figure 4.7: Weyl ratio for the level 5 Neumann eigenvalues. 
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function for both boundary conditions). Our data is inconclusive on this con- 
jecture. 

We are also interested in the difference between the two counting functions. 
By general principles, the Neumann counting function N^ N \t) will be greater 
than the Dirichlet counting functions N( D \t), and we expect 

NW(t)-N( Dy >(t)~t p (4.6) 

for some j3 that measures the relative sizes of 9SC and SC. A simple guess is 

/3 = ~ 0.4769. (4.7) 
logp 



In Figure |4.8| we graph the difference and in Figure |4.9| the associated Weyl 
ratio 

Again this is presumed to be asymptotic to a multiplicatively periodic function. 

In our data there is no evidence for large gaps in the spectrum, but there is 
evidence for spectral clustering (distinct eigenvalues that are close together, see 
[CSW ). The spectral asymptotics predicts that the separation between distinct 
eigenvalues should go to infinity on average. Nevertheless, in the Dirichlet spec- 
trum we have eigenvalue pair #53-54 with value 860.58 and #55 with values 
861.92, we have four eigenvalues #118-121 in the range 2005.8 to 2015.6, and 
four eigenvalues #152-155 in the range 2786.7 to 2787.8. 

Using our eigenfunction computations we were able to construct Dirichlet 
kernels 

N 

D N (x,y) = ^>2<Pj(x)<Pj(y) (4.9) 
i=i 

where {tpj} are the orthonormal eigenfunctions. To the extent that Djsr(-,y) 
resembles an approximate identity about the point this will tell us about the 
convergence of eigenfunction expansions. It was observed in |OSS| for the SG 
fractal that choosing TV so that the sum goes up to a spectral gap produces 
a strong approximate identity appearance, and later in [S2] it was shown for 
SG and related fractals that convergence of eigenfunction expansions for these 
special partial sums is much better than for ordinary Fourier series (continuous 
functions have uniformly convergent special partial sums). As already men- 
tioned, the SC spectra do not have such large spectral gaps, but even if we 
choose N to reach one of the more moderate gaps, the Dirichlet kernels do not 



show any approximate identity appearance. In Figure 4.10 we show the graph 
of Dn(x,v) for x restricted to a line segment in SC passing through the point 
y. The slow decay away from x = y is reminiscent of the Dirichlet kernels for 
ordinary Fourier series. 

We have also examined eigenfunctions and eigenvalues for periodic type 
boundary conditions corresponding to the construction of the torus, Klein bot- 
tle, and projective space from the square. SC with these identifications may be 
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N (N) (t>N (D) (t) 
beta- log 3/ log 1 0.01 1 =.4769 




Figure 4.8: The difference of the eigenfunction counting functions N^ N \t) 
N( D \t) at level 5. 



Neumann-DirichletWeyl Ratio (level 5) 
befca=log3Jlog1 0.01 1= 4769 




Figure 4.9: Weyl ratio for N^ N \t) - N^ D \t) given in p|. 
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Figure 4.10: Dirichlet kernel D^^x^y) on SC5 using the first 1200 Neumann 
eigenvalues where y is the cell '1555' (the center cell directly above the level- 1 
hole) and x varies along the horizontal line containing y. 



regarded as fractafolds based on SC in the spirit of [SI . Note that in the case 
of the projective space we create two singularities at the corner points of the 
square, identified in pairs. We show some of the results in Figures |4T[p~T3l and 
Tables E3E11 

There are some interesting coincidences in the Neumann spectrum and these 
spectra. For simplicity we denote the four cases by N, T, KB, and PS. The full 
D4 symmetry group acts on N, T, and PS, but not KB, and indeed there ap- 
pear to be no eigenspaces of dimension greater than one for KB. The one- 
dimensional eigenspaces for N and PS are identical. Exactly half the one- 
dimensional eigenspaces for N appear in the T and KB spectra, namely those 
that are symmetric under horizontal and vertical reflections (labels 1 + + and 
1 — h in Table 4.2). The two-dimensional eigenspaces for T and PS contain a 
one-dimensional subspace that contributes to the KB spectrum. (Note that the 
T and KB spectra also contain other one-dimensional eigenspaces unrelated to 
the other spectra.) 

All these coincidences are easily explained by symmetry considerations. In 
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Figure 4.11: The first nine nonconstant eigenfunctions on SC5 with torus bound- 
ary conditions 



A; 


Eigenvalue 


level m = 4 


level m — 5 


1 





2.10360108348969e-15 


-2.14570386198176e-17 


2 


32.2114556413112 


0.00320568102788810 


0.000320348757759280 


3 


32.2117876737192 


0.00320568102788826 


0.000320352059881687 


4 


33.0335325255338 


0.00328561228543238 


0.000328524461197705 


5 


46.1378042666663 


0.00458982380491410 


0.000458848815997367 


6 


67.3159845178793 


0.00670022323157356 


0.000669469652591202 


7 


70.3055712815539 


0.00699438207618352 


0.000699201634176303 


8 


70.3059096613708 


0.00699438207618378 


0.000699204999424840 


9 


91.1670279725552 


0.00906849340248511 


0.000906672597625728 


10 


112.219464252958 


0.0111377385936996 


0.00111604288766578 


11 


153.725854579008 


0.0152688242028740 


0.00152883145357489 


12 


153.726931440876 


0.0152688242028747 


0.00152884216316112 


13 


159.788304787400 


0.0158639440470136 


0.00158912361841407 


14 


164.837696052388 


0.0163805216075794 


0.0016393407286617!) 


15 


182.394416891886 


0.0181239623631475 


0.00181394549579459 



Table 4.3: The first 15 eigenvalues for levels m = 4 and 5 with torus boundary 
conditions 
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Figure 4.12: The first nine nonconstant eigenfunctions on SC5 with Klein bottle 
boundary conditions 





Eigenvalue 


level m = 4 


level m — 5 


1 





-3.02058725254462e-14 


-2.19737287176415e-17 


2 


32.2117498844683 


0.00320568102788416 


0.000320351684060745 


3 


33.032810672492 


0.00328561228541770 


0.000328517282238523 


4 


39.8922202086433 


0.00397107095374325 


0.000396735351869938 


5 


43.2714800226914 


0.00430597895848081 


0.000430342702485531 


6 


46.1377316931061 


0.00458982380490713 


0.000458848094240172 


7 


70.3057087137930 


0.00699438207617901 


0.000699203000964800 


8 


91.1664483134311 


0.00906849340247104 


0.000906666832810801 


9 


103.298679176177 


0.0102746017084947 


0.00102732406510132 


10 


112.219004510499 


0.0111377385936637 


0.0011160383154437 ( ) 


11 


122.693174548126 


0.0121921338975400 


0.00122020583266118 


12 


153.726542595912 


0.0152688242028903 


0.00152883829602755 


13 


162.162017612389 


0.0161071330316981 


0.00161273062218409 


14 


164.841195409480 


0.0163805216076038 


0.00163937553040158 


15 


182.388331131098 


0.0181239623631782 


0.00181388497180182 



Table 4.4: The first 15 eigenvalues for levels m = 4 and 5 with Klein bottle 
boundary conditions 
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Figure 4.13: The first nine nonconstant eigenfunctions on SC5 with projective 
plane boundary conditions 



A; 


Eigenvalue 


level m = 4 


level m — 5 


1 





2.15702836851335e-15 


-2.13740293909888e-17 


2 


17.7803169580100 


0.00177037771649474 


0.000176828471010163 


3 


33.0329798073433 


0.00328561228543052 


0.000328518964315239 


4 


39.8922635742169 


0.00397107095374375 


0.000396735783148419 


5 


39.8924984277846 


0.00397107095374659 


0.000396738118809666 


6 


46.1377416688106 


0.00458982380491383 


0.000458848193450360 


7 


72.9450515387370 


0.00724990726688279 


0.000725451743172754 


8 


76.4963443558820 


0.00760210249889050 


0.000760769993148170 


9 


91.1663022922193 


0.00906849340248375 


0.000906665380603414 


10 


112.220120732730 


0.0111377385937001 


0.00111604941647595 


11 


122.691748138782 


0.0121921338975364 


0.00122019164676202 


12 


122.692643256546 


0.0121921338975388 


0.00122020054887023 


13 


164.843452415657 


0.0163805216075811 


0.00163939797673661 


14 


178.014270428489 


0.0177288779201235 


0.00177038414625554 


15 


182.393894269051 


0.0181239623631491 


0.00181394029821588 



Table 4.5: The first 15 eigenvalues for levels m = 4 and 5 with projective plane 
boundary conditions 
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all the examples we may extend functions to be periodic of period 1 in x and y. 
In addition we have the following identities: 



• Case N: 

(Nl) u(-x,y) = u{x,y) 
(N2) u(x,-y) = u(x,y) 

• Case T: 

(Tl) u(x + l,y)=u(x,y) 
(T2) u(x,y + l)=u(x,y) 

• Case KB: 

(KB1) u(x + 1,1-7/) = u(x, y) 
(KB2) u(x,y + l) = u(x,y) 

• Case PS: 

(PS1) u(x + 1,1 -y)=u(x,y) 
(PS2) u(l-x,y + l)=u{x,y) 

In addition, if u is symmetric with respect to horizontal and vertical reflection 
then 

u(l — x,y) = u(x, y) and u(x, 1 — y) = u(x, y). (4.10) 



It is straightforward to show that under (4.10) all four sets of conditions are 
equivalent. This explains why the (1 + +) and (1 — h) N-eigenfunctions are also 
eigenfunctions on T, KB, and PS. On the other hand, if u is skew-symmetric 
with respect to horizontal and vertical reflections then in place of ( |4.10 ) we have 

— u(l — x,y) = u(x,y) and — u(x,l — y) = u(x,y). (4.11) 

Under this assumption we have 

u(x + 1,1 — 2/) = —u(—x, 1 — y) = u(—x, y) 

so (Nl) is equivalent to (PS1), and similarly (N2) is equivalent to (PS2). 

Thus the type (1 H — ) and (1 ) eigenfunctions in the N and PS examples 

are the same. 

Next consider a two-dimensional eigenspace for example T. In this eigenspace 
we can find an eigenfunction symmetric with respect to vertical reflections, 
u(x,l — y) = u(x,y). This makes (Tl) and (KB1) equivalent, and of course 
(T2) and (KB2) are identical equations, so this eigenfunction is also an eigen- 
function for KB. Similarly, in a two-dimensional eigenspace for PS there is an 
eigenfunction that is symmetric with respect to horizontal reflection, u(l — 
x,y) = u(x,y). This makes (PS2) and (KB2) equivalent, and (PS1) and 
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(KB1) are identical equations, so this eigenfunction is also an eigenfunction for 
KB. 

Of course the identical remarks are valid in the non-fractal case of the square 
with Neumann boundary conditions, the torus, the Klein bottle, and the pro- 
jective space, where they are presumably well-known, and can also be explained 
by the expressions for the eigenfunctions in terms of sines and cosines. 

Another interesting observation is the behavior of the differences between 
the eigenvalue counting function for these examples. In Figure 4.14| we show 
the difference between the eigenvalue counting functions for Neumann boundary 
conditions and the torus, showing a growth rate of t@ with /3 ~ 0.4769, and in 
Figure [47L5| we show the Weyl ratio. The se re sults are consistent with the results 
for 7V(^T(F) - N^ D \t) shown in Fi gures 4.8 and 4.9[similar results, not shown 



here, are obtained for the Klein bottle and projective space). In Figures 4.16 



and |4.17| we show the differences between the eigenvalue counting functions for 
torus and Klein bottle (respectively projective space). Here we see a dramati- 
cally different behavior, leading us to conjecture that the difference is uniformly 
bounded, (Note that we do not have enough data to distinguish between bound- 
edness and logt growth). Moreover there is rapid oscillation between positive 
and negative values. 

We note that for the standard Laplacian on the torus and Klein bottle, the 
difference of the eigenvalue counting functions is expected to be on the order of 



5 Heat kenrel 

The heat kernel h t (x,y) (for either Dirichlet or Neumann boundary conditions) 
is defined to be 

00 

h t (x,y) = ^2e- x i t (p j (x)(p j (y) (5.1) 
i=i 

where {(fj} is an orthonormal basis of eigenfunctions with corresponding eigen- 
values {Aj}. By integrating against the heat kernel 

u(x,t)= h t (x,y)f(y)dfjb(y) (5.2) 
J sc 



you can solve the heat equation 



— Q t U 0\ 

u(x,0) = f(x) ^ 



with the appropriate boundary conditions. Moreover, h t (x,y) gives the tran- 
sition probabilities of a Brownian motion path moving from x to y in time t. 
Using the probabilistic interpretation, Barlow and Bass [BB1993 were able to 
prove subGaussian bounds (above and below) for the heat kernel for small times, 

M.,»)-rW-3^!I) o«<i, («) 
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Neumann-Torus counting function. p~0.4840 




Figure 4.14: The difference between the eigenvalue counting functions for Neu- 
mann and torus boundary conditions. We estimate the exponential growth rate 
(shown with a red line) and obtain f3 « 0.4769. 



Neumann-Torus Weyl Ratio. p~0.4769 




Figure 4.15: The Weyl ratio for the difference between the eigenvalue counting 
functions for Neumann and torus boundary conditions with f3 « 0.4769. 
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Torus-Klein counting function. 




Figure 4.16: The difference between the eigenvalue counting functions for torus 
and Klein bottle boundary conditions. 



Torus-Projective counting function. 




Figure 4.17: The difference between the eigenvalue counting functions for torus 
and projective space boundary conditions. 
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where a = as before, and 7 and S are other constants. But this is not the 

log p ' ' 

end of the story. We would like to use our numerical data to shed some light 
on the more precise behavior of the heat kernel. In particular, the trace of the 
heat kernel (also called the partition function) 

00 

Z(t) = 5>- A '*. (5.5) 
i=i 

Figur e |5.1 1 shows the trace of the heat kernel with Neumann eigenvalues. Note 
that ( |5.4| ) implies that Z(t) ~ t~ a as t — >• 0. We obtain a slope of —a « —0.895 
which is reasonably close to our predicted a = ~ 0.9026. But again more 
is known, t a Z(t) converges to a multiplicatively periodic function bounded and 
bounded away from zero of period p as t — »• 0. In analogy with the case of SG 
we expect that t a Z(t) as a function of \ogt to resemble a sine curve with very 
small error. Figure [5T2] shows t a Z(t). We might also expect t a h t (x,x) for any 
point x to have a similar behavior, but the sine curve should show a phase shift 



that depends on x. Figure 5.3 shows t a h t (x,x) for four choices of x. 

Just as we studied the counting function and Weyl ratios of the difference 
between Neumann and Dirichlet boundary conditions, we do the same for the 
heat kernel. Again we expect Z^ N \t) — Z^ D \t) ~ where /3 = log 3/ log p « 
0.4769. Figure [5~4] shows the graph of the difference of the heat kernels and has 



an approximate slope -0.4797 which is satisfyingly close to —f3. If we exa mine 



t^{Z^ N \t) — Z( D \t)) we see a similar periodic behavior as shown in Figure 5.5 



Figures [5^6| ?? shows graphs of the off-diagonal heat kernel ht(-, y) for a few 
choices of fixed y. Note in particular that the level sets in Figure ?? resemble 
the level sets of the resistance metric in Figure [3T8| even though the functions 
are so different. 



6 Boundary Behavior 

It is easy to obtain the analog of the Gauss- Green formula for functions u 
satisfying Neumann boundary conditions. We simply write 

S m (u,v) = -^^2v(x)^2(u(x)-u(y)) (6.1) 

x y<^x 

= -E 8 -m*)EG 

x yr^x ^ 

= -^8-^(x)A m ^(x). 



(u(y) - u{x)) 



If u satisfies Neumann boundary conditions then we obtain 

£(u,v) = — / vAud/i (6.2) 

JSC 
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(a) 



(b) 



10' 



to: 




10 1 10 2 10 3 10" 



(e) 

Figure 5.3: t a h(x,x) on a log-log scale using the first 1500 level-5 Neumann 
eigenvalues and eigenfunctions. In (a) x = '11 111 7 , the top-center cell; (b) 
x = '00000', the top- left corner cell; (c) x = 15555, the center cell immediately 
above largest hole of SC; (d) x = '10GSB', the level- 1 junction point along the 
top edge of SC; (e) x = '05271 7 , an arbitrary interior cell. 




Figure 5.4: Difference of Neumann and Dirichlet traces of the heat kernel 
Z^ N \t) — Z^ D \t) on a log-log scale using all level-5 eigenvalues. The slope 
for small t is -0.4797 (shown in red). The apparent change of behavior for very 
small values of t is a simple consequence of the fact that we have not computed 
enough eigenvalues. 




Figure 5.5: t^{Z^ N \t) — Z^ D \t)) on a log-log scale using all level-5 eigenvalues. 
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Figure 5.6: Graphs of the off-diagonal heat kernel h t (-,y) for y fixed as the 
corner point '00000' with t = 1/100 (left) and t = 1/10 (right). The first row 
shows h t (-,y) in its entirety. The second row shows the restriction of h t (-,y) to 
the horizontal line containing y. The third row shows the kernel with level sets 
colored differently. 
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Figure 5.7: Graphs of the off-diagonal heat kernel h t (-,y) for y fixed as the 
center point '15555' and t = 1/100 (left) and t = 1/10 (right). The first row 
shows h t (-,y) in its entirety. The second row shows the restriction of h t (-,y) to 
the horizontal line containing y. The third row shows the kernel with level sets 
colored differently. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 5.8: Graphs of the off-diagonal heat kernel h t (-,y) for y fixed as the top 
point '11111' and t = l/100(left) and t = 1/10 (right). The first row shows 
h t (-,y) in its entirety. The second row shows the restriction of h t (-,y) to the 
horizontal line containing y. The third row shows the kernel with level sets 
colored differently. 
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in the limit as m — >• oo. If u does not satisfy Neumann boundary conditions 



then we do not have the correct expression for A m u(x) in (6.1) when x is a 
boundary cell, because we need to take into account the values when y = x*, 
the virtual neighbor (two of these at corner points). 
Thus we have 

f m (ii^) = -^8- m ^(x)A m ^(x)+ v{x)p- rn {u{x)-u{x^). (6.1') 



X xed SCrn 

In the limit we hope to obtain 

£(u, v) — — I vAud/i+ / vd n udu (6-2') 

JSC Jd SC 

for some appropriate measure on the boundary, and some appropriate definition 
of the normal derivative d n u. Since the boundary of SC is the same as the 
boundary of the square, one might speculate that dv could just be Lebesgue 
measure, but our data does not support this supposition, and indeed there is 
no reason to believe that the measure should be so place independent, since the 
geometry of SC near boundary points varies considerably from point to point. 

We are not prepared to put forward a conjectural definition of the normal 
derivative. Instead we investigate the preliminary question: what is the behavior 
of u(x) as x approaches a boundary point zl When u(z) = 0, as it does in a 
number of our computed examples, this becomes the question of the decay rate 
of u{x) as x — ^ z. 

First, we will outline the methods used to measure decay rates. Let z be a 
point on <9SC n . Suppose there exist three n cells in a line perpendicular to the 
edge containing z. We denote them by Xj for j = 1,2,3 with z G x\. Thus, 
we exclude points z that lie near an n-level hole and for any boundary edge 
of SC n , we can calculate decay rates for exactly |3 n of the n-cells intersecting 
that edge. Note that the average distance to the boundary for the cell Xj is 
^~y 2 ^ for j = 1,2,3. To find the decay rates we will be measuring at what 
rate u(xj) — » for j = 1, 2, 3, when u{z) = 0. 

A naive guess would be that 

u(x j )-u(z)^A ^ J ~ 3 l /2) y (6.3) 

for some power a and A ^ 0. In that case, we would set A = d n u(z). 

We present decay rates of various functions that vanish along <9SC starting 
with the harmonic functions from Figure 3.1 , namely, the harmonic functions, 
hk, defined by sin(fc7r£) along the top edge and zero along the other edges along 
9SC. We will look at the decay rates given by ( |6.3| ) where z is along the bottom 
edge (furthest from sin(fe7rt)). We use the first 6 harmonic functions, hk, at the 
highest level available (for us, level-6 harmonic functions). For each of the | • 3 6 
acceptable cells along the bottom boundary, we can compute estimates for A 



and a in (6.3). 
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Figure 6.1 



shows estimates for a for each of the | • 3 6 cells along the bottom 
boundary of SC6- The graphs of a for hk for k odd are nearly identical to each 
other; the same is true for k even (they differ by at most 10 -4 ). It is clear from 
the data that a is not constant, but depends on the point z, with an average 
value of 1.1135. The behavior for z near a corner point is different as might be 
expected, but otherwise the variation of a with z follows a distinctive pattern. 

Figure [6^2] shows the estimates for A for each of the hf.. Just as with our 
estimates for a, the graph shapes for k odd and even are virtually indistinguish- 
able. 

We next present the same results for Dirichlet eigenfunctions. We present 
data using the first six level-5 Dirichlet eigenfunctions ipk- We compute a and 



A exactly as was done for the harmonic functions above. Figure 6.3 shows the 
estimates for a for the first six cp^. Just as with the harmonic functions, the 
behavior of a with z is similar with the exception of the corners and the centers 
for the even ip^. These estimates for a vary slightly more than for the harmonic 
functions but we obtain an average value of a w 1.1152. This is slightly higher 
than the a ~ 1.1135 but this is most likely due to the higher variation at the 
corners. Figure 6.4 shows the estimates for A. Just as with the harmonic 
functions, the general shape for the (p dd and (p even are preserved. Figure |6.5 
shows the simultaneous graphs of a for the four examples hi,h2,<pi,y>2- 

As we saw, the boundary behavior acts quite differently at the corners of 
<9SC. We now focus on corner behavior and present some results. 

We calculate the convergence at corner points exactly as before except for 
j = 2 and 3, u(xj) becomes the average value of the two cells that are (j — 1) 
cells away from z C x\. That is, we calculate the corner decay rates from the 
two perpendicular directions at the corner. Again, we guess that 



u(xj) — u(z) « A 



(i - 1/2) 



First, take the harmonic functions from before, h^. There are two types 
of corner points for these functions. There are 'top' corners, which have one 
edge equal to zero and one with sin(7rA:x); and there are bottom corners, which 
have both edges set to zero. For this reason, it is natural to expect the bottom 
corners to exhibit a higher value of a which we observe as shown in Tables 6.1 
and 16.21 Also note that each of the harmonic functions have the same 'bottom' 





k 




1 


2 


3 


4 


5 


6 


a 


0.6884 


0.6882 


0.6879 


0.6875 


0.6871 


0.6867 


A 


0.2802 


0.5595 


0.8374 


1.1134 


1.3872 


1.6587 



Table 6.1: Values of a and A for the top corner of the harmonic functions hk. 
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(e) h 5 42 
Figure 6.1: Estimates of a for h^. 





(a) hi 



(b) h 2 




0.4 0.6 0.8 




(c) h 3 



(d) h 4 




0.4 0.6 0.8 



(e) h 5 
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(f) h 6 



Figure 6.2: Estimates of A for h^. 





0.2 0.4 



(a) ip! 



(b) (p 2 





0.2 0.4 



(c) ^3 



(d) <p 4 





0.2 0.4 



(e) ^5 
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(f) ^6 



Figure 6.3: Estimates of a for ip^. 




0.2 0.4 0.6 0.8 




0.2 0.4 0.6 



(a) LP! 



(b) (p 2 





0.2 0.4 0.6 



(C) ^3 



(d) <p 4 





0.2 0.4 0.6 



(e) ^5 
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(f) ^6 



Figure 6.4: Estimates of A for 




1 

0.95 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




Figure 6.5: Estimates for a for hi,h,2,ipi, and cp2 for the whole interval and 
zoomed in to a small interval away from the corners and center point. Note 
that in the zoom the four graphs are virtually indistinguishable. 
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k 




1 


2 


3 


4 


5 


6 


a 


2.6903 


2.6903 


2.6903 


2.6903 


2.6903 


2.6903 


A x 1(T 3 


0.1356 


0.0623 


0.0077 


0.0178 


0.0207 


0.0005 



Table 6.2: Values of a and A for the bottom corner of the harmonic functions 
hk- 

We do the same computation for the Dirichlet eigenfunctions. However, for 
these functions, we do not need to test top and bottom corners. Eigenfunctions 
corresponding to an eigenvalue with multiplicity 1 will automatically be sym- 
metric from any corner. Eigenfunctions corresponding to an eigenvalue with 
multiplicity 2 will have both corners accounted for since the two eigenfunctions 
are just a 90 degree rotation of one another. All values in Table [63] are taken 
from the top left corner of SC. 





k 




1 


2 


3 


4 


5 


6 


a 


1.4119 


1.4119 


1.4119 


1.4119 


1.4119 


1.4119 


A x 10- 4 


0.3131 


0.3240 


0.2795 


0.1298 


0.1881 


0.1385 



Table 6.3: Values of a and A for the corner of the Dirichlet eigenfunctions (p^. 



7 Periodic covering spaces 

There are many covering spaces of SC. We are interested in spaces which have 
a Z-action of isometries, for then we may use ordinary Fourier series to describe 
the spectrum of the Laplacian in terms of eigenfunctions on SC satisfying certain 
periodic boundary conditions. Analogous examples for Julia sets are discussed 
in [ADS]. We discuss two examples. The first is the staircase in which the eight 
1-cells of SC are cut along one edge and the cut edges are joined up to the next 
level (up and down) of a spiral staircase. The Z action moves up and down 
the levels, and the projection onto SC just maps each 1-cell of the staircase to 
the corresponding 1-cell in SC. The second example, the strip, consists of an 
infinite horizontal row of copies of SC identified along their vertical boundary 
edges. This is strictly speaking not a covering space because there is no covering 
map, but it is close enough in spirit. It is a fractafold [Si], in that every point 
has a neighborhood that is isometric to an open set in SC so in particular the 
definition of the Laplacian lifts to the strip. The Z-action is just horizontal 
translation. 

We describe the situation for the strip first, since it is simpler. The strip is 
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a subset of R x [0, 1]. For any in [0, 1] we look for eigenfunctions satisfying 



u(x + l,y) = e'™°u(x,y) 



(7.1) 



and Neumann boundary conditions on the top and bottom lines of the strip. 
We call this the ^-spectrum. To be specific, on level 1 we have the eight 1-cells 



xq,...,x? (see Figure 7. 



x 


Xi 


x 2 


x 7 




x 3 


X 6 


X 5 


£4 



Figure 7.1: The eight level- 1 cells of SC used in describing the strip covering 
space. 

At x\ and x 5 the eigenvalue equations are the same as for SC 

2u(x\) — u(xq) — u(x 2 ) — \u[x\) 
2u(x$) — u(x±) — u(xq) = Xu(x^) 

However, the other 1-cells have neighbors in the strip that are translates by 
±1 in the x direction of other 1-cells, so 



3u(x ) 


- u( Xl ) - u(x 7 ) - e~ 27Tii 


} u(x 2 ) 


— Xu(xq) 


3u(x 2 ) 


- u( Xl ) - u(x 3 ) - e~ 27rU 


] u(x ) 


= Xu(x 2 ) 


3u(x 3 ) 


- u(x 2 ) - u(x 4 ) - e~ 27TU 


} u(x 7 ) 


\u(x 3 ) 








etc. 



Note that the complex conjugate of a #-eigenfunction will be a (1— #)-eigenfunction, 
so the ^-spectrum is real and invariant under reflection about = \. Vertical 
reflections about y = \ in the strip maps ^-eigenfunctions to ^-eigenfunctions, 
so the ^-spectrum splits into the even and odd eigenfunctions. Horizontal re- 
flections map ^-eigenfunctions to (1 — #)-eigenfunctions. Since the symmetry 
group of the strip is abelian, we do not expect multiplicities greater than one in 
any of the ^-spectra, except for coincidences. 

The spectrum of the Laplacian on the strip is continuous, and consists of 
the union of all the 0-spectra. In Figures f7T2]|7.3| we show graphs of the bottom 
of the ^-spectra, the horizontal axis being and the vertical axis being A, the 
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^-eigenvalues. It is easy to identify the vertical symmetry type of the eigen- 
functions for = and = \. Of course the ^-spectrum varies continuously 
with so we see a set of curves going from 6 = to = 1 with some crossings 
and near crossings. Generically, we do not expect to see crossings for pairs of 
eigenfunctions of the same symmetry type, but crossings may occur for pairs 
with opposite symmetry types. It is difficult to distinguish crossings and near 
crossings from the data, but tracking symmetry types is helpful in deciding. See 



Figure 7.4 



In Figure |7.5| we show the bottom portion of the spectrum of the strip, 
obtained by projection the ^-spectra onto the A-axis (shown horizontally in the 
figure). We observe that it is a union of intervals. It is not clear whether or not 
we should expect this pattern to persist for larger A values, or if the spectrum 
contains an infinite interval [a, oo). 

For the staircase, we assume that the cut is placed between 1-cells xj and 
xq. The ^-eigenvalue equations on level 1 are the usual away from the cut, 

2u(xj) — u(xj-i) — u(xj+i) — \u(xj), j = 1, ...,6, 

and 

2u(x ) - u(xi) - e~ 27ri9 u(x 7 ) = Xu(x ) 
2u(x 7 ) - u(x 6 ) - e 27Tl9 u(x ) = Xu(x 7 ) 



along the cut. Again the complex conjugate interchanges ^-eigenfunctions and 
(1 — #)-eigenfunctions, so we have ^-spectra symmetric about 6 — |. 

In this example we have an additional symmetry of \ rotation, denoted R. 
(Note that we could have used R to generate a more refined Z-action, but this 
would require a choice of fundamental domain that does not split into a simple 
union of 1-cells.) For a #-eigenfunction we have 

u(R A x) = e 2 ^ e u(x) (7.2) 



as the analog of (7.1), so we can further split the eigenspaces according to the 
symmetry types 

u(Rx) = e 27Ti ^^u(x), k = 0, 1, 2, 3. (7.3) 

For generic there will be four independent families (for each choice of fc), 
but for 6 — or — \ there will be families of multiplicity two. Of course = 
corresponds with the Neumann spectrum of SC, and the cases 

u(Rx) = u(x), for/c = 0, (7.4) 

and 

u(Rx) = —u(x), for k = 2, (7.5) 



49 



0.25 




I i i i i_ 

0.2 0.4 0.6 0.8 



Figure 7.2: Graph of the lowest twenty eigenvalues of the ^-spectrum for the 
strip covering space with the horizontal axis being 9 and the vertical axis being 
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0.2 



0.4 



0.6 



0.8 



Figure 7.3: Graph of the lowest twenty eigenvalues of the ^-spectrum for the 
staircase covering space with the horizontal axis being 9 and the vertical axis 
being A. 
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Figure 7.4: First 12 eigenvalues of the spectrum of the strip with the symmetry 
types labeled. 
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0.2 0.4 0.6 0.S 1 1.2 1.4 1.6 



Figure 7.5: The ^-spectra of the strip projected onto the A axis. It is a union 
of intervals. 



give the multiplicity one spaces, while the k = 1,3 cases combine to give a 
multiplicity two space satisfying 



u(Rx) = —v(x) 
v(Rx) = u(x) 



(7.6) 

For 6 = \ we can group the four cases into groups of two as follows: 



u(Rx) = e™ u(x) = ^(1 + i)u(x), k = , , 

u{Rx) = eT U (x) = |(l-i)t/(4 k = 3 



{u(Rx) = e 3 ^u(x) = ^(-1 + i)u(x), k = l , 
u(Rx) = eT^) = |(-l-i)u(x), k = 2 ' 

Replacing the complex valued u by u + iv for real- valued u and v yields 



I u(Rx) = ^(u(x)-v(x)) 

and 



2 »y-JJ ( 7J> ) 



f u(ite) = &{-u{x) - v{x)) , . 

\ v(Rx) = ^(u(x)-v(x)) ' 1 ' ^ 

Thus at = \ all eigenspaces have multiplicity 2 (or a multiple of 2), and 
they come in two symmetry types that we denote (0,3) and (1,2). Moreover, 



as varies between and |, the multiplicity two eigenspace of type (7.6) will 
split and one curve will join up at = \ with one curve starting from ( |7.4 ) to 



form a (0,3) eigenspace, while the other curve will join up with one from (7.5) 
to form a (1,2) eigenspace. Thus we expect to see groups of 4 curves (or 4n if 
there are crossings) in the ^-spectra. This is apparent in Figure 7.3 , where we 
see groups of 4,4,4,4,8,... starting from the bottom, and sizable gaps between 
groups. Figure |7.6| shows the spectrum of the staircase, and in this example it 
appears likely that the pattern will persist and yield an infinite disjoint union 
of finite intervals. 
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1 2 3 4 5 

Figure 7.6: The ^-spectra of the staircase projected onto the A axis. It is a 
union of intervals. 

8 Conclusions 

In this paper we have presented a substantial body of numerical computations 
of basic functions on SC. While there is no theoretical estimate for the rate of 
convergence of the algorithms we used, the results obtained give strong exper- 
imental evidence that we are in fact obtaining reliable accuracy. For some of 
the computations it would have been desirable to continue to greater depth in 
the approximation to SC, for the most part the error is acceptably small. So we 
achieved a "picture book" glimpse of basic calculus on SC. 

One general conclusion is that calculus on SC is somewhat more regular 
than on SG and other PCF fractals. Functions such as harmonic functions, 
eigenfunctions of the Laplacian, Poisson kernels, and heat kernels appear more 
regular when graphed against the Euclidean embedding of SC in the plane. 
The spectrum of the Laplacian, as observed in the Weyl ratio of the eigenvalue 
counting function, is somewhat smoother. L 2 -normalized eigenfunctions of the 
Laplacian appear to remain uniformly bounded throughout the spectrum, or at 
worst have logarithmic growth with the eigenvalue. On the other hand, since 
there are apparently no large gaps in the spectrum, there are no nice Dirichlet 
kernels with approximate identity behavior, so partial sums of Fourier series are 
not well-behaved. 

The geometry of SC is very far from being homogeneous. In particular there 
are "corner points", lying at the corners of each of the deleted squares in the 
construction of SC. These points have the same local geometry when the scaled is 
reduced, as compared with generic points, where the local geometry varies with 
the scale. It is not surprising that the behavior of functions at corner points is 
somewhat different from the behavior at generic points. Further work is needed 
to elucidate these differences. But perhaps the difference is not as pronounced 
as in the junction point /nonj unction point dichotomy for PCF fractals. 

What is the "natural" metric on SC? There are several competing metrics 
to consider: the metric inherited from the embedding of SC in the plane, the 
geodesic metric inherited from that embedding, and the effective resistance met- 
ric. There are also more qualitative geometric objects, like the level sets of the 
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heat kernel, that are roughly comparable to spheres in the effective resistance 
metric. We would need to be able to compute effective resistances to greater 
resolution in order to draw more precise conclusions. 

We have also been successful in understanding the spectrum of several fractafolds 
based on SC. In the compact case we have looked at boundary identifications 
corresponding to those that produce the torus, Klein bottle, and projective 
plane when applied to the square that contains SC. Our most striking conjec- 
ture based on experimental evidence is that the differences of the eigenvalue 
counting functions for these three fractafolds remains uniformly bounded, a 
statement that is evidently false in the nonfractal case. In the noncompact 
case we have looked at two periodic covering spaces, the strip (not technically 
a true covering space) and the staircase. We have experimental evidence that 
the spectrum is continuous and consists of a countable union of intervals. 

We have not been successful in defining a normal derivative on the boundary 
of SC with a corresponding Gauss- Green formula. However, we have presented 
a wealth of data on the decay rates of functions near the boundary, and this 
should be useful in further work on this interesting and important problem. 
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